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ABSTRACT 


We  derive  the  closed  form  expression  for  the  bit  error  probability  of  dense  WDM 
systems  employing  an  external  OOK  modulator.  Our  model  is  based  upon  a  close 
approximation  of  the  optical  Fabry>Perot  filter  in  the  receiver  as  a  single-pole  RC  filter  for 
signals  that  are  bandlimited  to  a  fiequency  band  approximatdy  equal  to  one  sixtieth  of  the 
Fabry-Perot  filter's  free  spectral  range  Our  model  can  handle  bit  rates  up  to  2.5  Gb/s  for 
a  free  spectral  range  of  3800  GHz  and  up  to  S  Gb/s  when  the  power  penalty  is  1  dB  or 


less. 
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L  INTRODUCTION 


Wavelength  division  multiplexing  (WDM)  systems  have  been  increasingly  proposed 
as  an  attractive  alternative  to  coherent  optical  frequency  division  multiplexing  (FDM) 
systems  [Ref  l-S],  Although  WDM  systems  with  direct  detection  do  not  have  the 
channd  aq)acity  of  coherent  optical  FDM  systems,  they  are  much  less  costly  to 
implement.  Furthermore,  present  fiher  technology  enables  the  designers  to  tightly  pack 
the  channel,  resulting  in  dense  WDM  systems  that  can  provide  aggregate  bit  rates  of  numy 
terabits  per  second  (1  Tb/s  =10'*  b/s).  Dense  WDM  systons  are  particularly  attractive  in 
the  area  of  undersea  surveillance  where  hundreds  of  sensors  and  data  collection  sites  are 
envisioned  being  merged  onto  single-fiber  superhighways  through  massive  data  fusion 
Other  applications  call  for  relatively  low-bandwidth  data  collection  over  many  months  to 
be  dumped  quickly  to  a  remote  recording  site  in  a  matter  of  minutes.  This 
"collection-and-dumped"  compression  can  demand  total  data  rates  on  the  order  of 
hundreds  of  Gb/s.  Long  distance  between  data  collection  sites  and  a  remote  recording  site 
requires  the  use  of  optical  amplifiers.  Therefore,  it  is  necessary  to  pack  all  channels  within 
the  optica]  amplifier  bandwidth. 

A  dense  WDM  receiver  with  on-ofif-keying  (OOK)  modulation  can  be  modeled  as 
shown  in  Fig.  1.  Conceptually,  the  analysis  involves  two  main  operations:  1)  a 
convolution  operation  to  evaluate  the  signal  at  the  output  of  the  optical  filter,  a 
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Fabry>Perot  (FP)  filter  in  our  investigation,  and  2)  the  integration  of  the  output  of  the 
photodetector  Evaluation  of  bit  error  probability  by  the  numerical  analysis  of  these  two 
operations  has  been  carried  out  in  [Ref  6],  with  a  number  of  approxinutions  nuuie  to 
reduce  the  computational  complexity.  In  this  investigation  the  FP  fiher  is  shown  to  be 
well  approximated  by  an  RC  filter  within  the  fi’equency  range  |/”/o  |  <  FSR  /  207t,  where 
FSR  is  the  free  spectral  range  of  the  FP  filter  [Ref  7]  and  is  the  FP  filter  center 
frequency  For  example,  given  FSR  =  3800  GHz,  the  approximation  works  very  well  for 
I  f -fo  I*'  GHz;  that  is,  the  effects  of  adjacent  channels  within  a  121  GHz  bandwidth 
centered  at  /,  must  be  included,  while  all  others  can  be  neglected.  This  simple  model 
agrees  well  with  [Ref  6]  as  demonstrated  in  Section  HI  Furthermore,  this  model  enables 
us  to  obtain  a  closed  form  analytical  expression  for  the  bit  error  probability  for  which 
numerical  results  can  be  obtained  with  little  effort.  Our  investigation  shows  that  this 
simple  model  provides  accurate  results  as  compared  to  those  in  [Ref  6]  for  bit  rates  up  to 
2.5  Gb/s  when  the  effects  of  four  adjacent  channels  are  included  with  FSR  =  3800  GHz 
Actually,  when  the  power  penalty  relative  to  single  channel  operation  is  1  dB  or  less,  there 
is  virtually  no  difference  in  the  effect  of  four  or  two  adjacent  channels  Thus,  for  this 
power  penalty  criterion,  this  simple  model  can  handle  bit  rates  up  to  5  Gb/s  for  a  FP  filter's 
FSR  =  3800  GHz. 

In  Section  II  the  closed  form  expression  for  the  decision  variable,  and  consequently, 
the  bit  error  probability  assuming  all  channels  are  bit  synchronous  as  in  [Ref  6]  is  derived 
Section  III  presents  the  numerical  results  which  include  the  bit  error  probability  versus  the 
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signal-to-noise  ratio  as  function  of  the  FP  filter  bandwidth  and  channel  spacing,  and  the 
power  penalty  (relative  to  single  channel  operation  without  filtering  or  with  filtering  but  no 
intersymbol  interference)  versus  the  channel  spacing  as  a  function  of  the  bandwidth 
Finally,  a  summary  of  results  appears  in  Section  IV 
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a  ANALYSIS 


The  receiver  model  for  the  dense  WDM  system  is  shown  in  Fig.  1.  The  desired 
signal  is  filtered  by  a  Fabiy-Perot  (FP)  filter  that  rqects  adjacent  channels  The 
photodetector  is  assumed  to  have  a  responsivity  £  (AAV).  The  detected  current  is 
amplified  by  a  low  noise  amplifier  that  contributes  a  postdetection  thermal  noise  n{t)  with 
spectral  density  Ng  (A‘/Hz)  The  decision  variable  at  the  output  of  the  integration  is 
compared  to  a  threshold  a  to  determine  whetha*  a  bit  zero  or  bit  one  was  present 
A.  INPUT  SIGNAL 

For  convenience,  we  designate  channel  0  as  the  desired  channel,  and  channel  k  as 
an  adjacent  channel  where  k  =  -A//2,  . ..,  -1,  1,  ...,  M2  with  M  an  even  integer  We 
consider  the  equivalent  lowpass  (complex  envelope)  data  signal  in  channel  0  and  channel  k 
as  follows: 

bo{t)  =  bojpTit-iT) 

i=-Lo 

bk(t)  =  bkjeJ^'^'pjit-lT) 

where 

T:  bit  duration 

K,  =  {0, 1 } :  bit  in  channel  0  in  the  time  interval  {iT,  (i+l )7) 

^k.1  ={0,e^*‘>is  the  A*  bit  in  channel  k  in  the  time  interval  (/7’,(/+l)7) 


(1) 

(2) 
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Figure  1  OOK  receiver  structure 


a  phase  ofl^  between  channel  k  and  channel  0  and  is  assumed  to  be  uniformly 
distributed  in  (0,  27c)  radians 

00,^  radian  frequency  spacing  between  channel  k  and  channel  0  with  0^  =  -<o^ 

The  function/?^/-//)  is  defined  as 


iT<t<{i+\)T 

otherwise 


(3) 


In  both  (1)  and  (2),  the  non-negative  integers  and  L  represent  the  number  of  bits  in 
channel  0  and  k,  respectively,  that  proceed  the  detected  bits  ^The  received  dense 
WDM  equivalent  lowpass  signal  at  the  input  of  the  FP  filter  is  given  by 

Hi)  =  JPbHty^  f  JPbHO  (4) 

k=-MI2 
k*0 


where  P  is  the  received  optical  power. 

B.  FABRY-PEROT  FILTERED  OUTPUT  SIGNAL 


The  FP  filter  can  be  characterized  by  the  following  equivalent  lowpass  transfer 
function  [Ref  1,7] 


//(/)  = 


1-p 

l_pg-j2,fFSR 


\-A-p 

1-p 
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(5) 


«(/)  =  I — —w, 

l-pcos(^)+;psin(^) 


l-A-p 
*  1-p 


where  p  is  the  power  reflectivity,  A  is  the  power  absorption  loss  (zero  for  an  ideal  FP 
filter)  and  FSR  is  the  free  spectral  range.  For  \f  <  FSR/lOit  and  assuming  i4  =  0,  we  can 
approximate  Hi/)  as  follow: 


w(/)» 


1-p 


P^^fSR  i\-tt)FSR 


,  w 

1+7-r 


,\/\<FSm0n 


(6a) 


where 


FSR{l-p) 

P 


(6b) 


The  fi^ee  spectral  range  FSR  can  be  related  to  the  full  width  at  half  maximum  (FWHM) 
bandwidth  B  and  the  finess  F  of  the  FP  filter  as 


FSR  = 


I-p 


=  5F 


(7) 


Thus  if  the  signal  is  bandlimited  to  |/|  <  FSR/20k,  we  can  truly  approximate  (5)  with  a 
single>pole  RC  filter  with  the  following  transfer  function  and  impulse  response 


//(/)  = 


(8) 
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h{t)  -  ce~^‘ ,  /  >  0 


(9) 


Figures  2a-b  show  the  magnitude  and  phase  (radians)  of  H{f)  of  the  FP  fiber  in  (S)  and 
its  single-pole  RC  filter  approximation  given  in  (8)  for  p  =  0.99,  F=3 12.6, 5=121  GHz 
and  FSR  =  3800  GHz.  Note  that  as  the  frequency  increases,  the  phases  of  the  FP  filter 
and  the  RC  filter  differ  markedly,  but  the  magnitudes  of  their  transfer  functions  ronain 
identical  and  attenuate  rapidly  When  |  /  |  >  FSR/lOn,  the  magnitude  of  Hi  f)  is  very 
small,  and  therefore,  the  effect  of  adjacem  channel  interference  beyond  this  frequency 
range  is  negligible.  Figure.  3  shows  the  normalized  impulse  response  of  both  FP  and 
single-pole  RC  filters  In  summary,  the  above  ^proximation  is  valid  for  dense  WDM 
analysis  when  the  filter  finess  F  is  large  or  equivalently  the  FWHM  bandwidth  B  is  small 
since  the  equivalent  lowpass  signal  must  be  bandJimited  to  about  |  /  |  <  FSR/lOn. 

This  approximation  has  been  used  in  [Ref  5]  to  study  spectral  efficiency  of  optical 
FDM/ASK  systems,  which  involves  the  evaluation  of  the  decision  variable  for  worst-case 
analysis  using  the  eye  diagram  technique.  Since  we  are  interested  in  the  detected  bit  b  g  ^ 
in  the  time  interval  (0,7'),  we  consider  the  output  filtered  signal  s(  t),  0<  t  <  T  given  by 

s{t)  -  Ssif)  +  Sisiif)  +  SACj{t)  ,  0  <  f  <  r  (10) 


where 

Sg  (r):  desired  signal 

s.^it ):  intersymbol  interference  signal 
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Phase  of  H(f)  ( 


(a) 


Figure  2:  Spectral  characteristics  of  the  Fabry-Po'ot  filter  and  the  «4}proxiniated 
single-pole  RC  lowpass  filter 
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time  (ns) 


Figure  3:  Normalized  impulse  response  of  the  Fabry-Perot  filter  and  the  iqiproximated 
single-pole  RC  lowpass  filter 
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(11) 


sb(0  =  JPbQ,Q  J  h(J-i)dx 
0 


=  o<f<r 


Sisiit)=jp  t  bo^i  J  hit-x)dx 

t=-Lo  iT 


=  JPe-“ 


e~"  z  bo,,(e<‘*'>‘’'-e‘‘’^),0<t<T 

i=—Lo 


(12) 


sacM^Jp  i  {[  i  **./  J  /i(r-T)e^‘'rft] 

k=-M/2  l=-L  If 


k=-M/2  l=-L 
k*0 


+i*,o  j  *('-  T)e^‘’rft},  0  <  /  <  r 

0 


JPce-"  T  rir{[  ^  6*., 


ik=-A<y2  /=-! 

/M) 


-e<‘*>‘)'^)]  +  *t,o(e<"'^‘)'-l)>,  0<r<7'  (13) 

The  FP  filtered  output  s{t)  is  detected  by  the  photodetector  which  produces  a 

current  of  /f  |  s(t)  p  An^s.  This  current  plus  additive  white  postdetection  thermal  noise 

current  from  the  amplifier  is  integrated  by  the  imegrator  to  obtain  a  decision  variable  for 

the  threshold  detector. 
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C.  DECISION  VARIABLES 

The  decision  variable  T  appearing  at  the  integrator  output  consists  of  the  figwal 
component  X  and  noise  component  V 

r  =  (14) 

where 

X=\  ^\s{t)\^dt  (15) 

0 

T 

n{t)dt  (16) 

0 

We  note  that  is  a  zero  mean  Gaussian  random  variable  with  variance  NJ.  Substituting 
(10)-(13)  into  (15)  we  obtain  the  signal  component  A' as  a  function  of  the  two  parameters 
cT  and  (nj,  which  represent  the  effect  of  intersymbol  interference  and  adjacent  channd 
interference,  respectively 


£PTbl,[\-jf{\  -e-0  +  s?(l 


2cT' 


t=-Lo 


k=-Mll  h-L  ^+J<i*kTtcT 


^g(cT+j&kT)(l+l)  _^(cr-f7o*r)/^|2 
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k=-MI2 

k*0 


f 

m=-Mn. 

nt^ 


jii,0 


( X-^-juikT !  cT)  ( \-jvimT  /  cT) 


ffCT,  (0fc=<0m 
I L 

cT  ’ 


l-jfojtT/cT 


l+JHimTIcT 


^S,I  ¥  ¥ 


A=nMy2  m=-M/2 
k*Q  m*0 


i 


_ **./*1.0 _ 

( l+yco*r  /  cTX  X-ji&mTIcT) 


4(1  -  1  >  +  -  2e-"0 


-2.-)  T 

k=-Mf2 

M 


**.; 

X+j&kTIcT 


^g{cT+jakT){M) 


_^(cr+7«)tn/)  +  2 

ifc=-Af/2 


**■0  /  .  --cf 

X+jfatT  IcT  ^jiaicTIcT 
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l-JtaitTIcT 


)}+^  i  bo,(e^‘^'^^ 

i=-Lo 


,icP 


)Re{ 


f  -i 


^k.l 


k=-MI2 

k*Q 


l-_l  \-kj<sikTlcT 


(1 


^^(cr+;fi)krx/+l)  _  gicT-¥jfSikT)l-^ 


k=-MI2 

M 


bk.0  2cr 

l+jOkTIcT^  X-jiSikTIcT 


(17) 


D.  BIT  ERROR  PROBABILITY 

For  a  detection  threshold  a  and  an  ISI/ACI  bit  pattern  l>~{b^,  },  i  =  -1; 

/  =-£,  .  ,  0;  ib*  44(2,  ...,  M/2,  k*0,t!t»  conditional  bit  error  probability  of  the  OOK 
signal  represented  by  the  Gaussian  random  variable  Y  in  (14H16)  is  given  by  [Ref  8] 


where  Q  (X)  is  defined  as 


(18) 


e(jt)  =  -^  J  e->''^dy  (19) 

and  A'o  and  A',  are  the  values  of  AT  in  (17)  for  =  0  and  bo,o  ~  1«  respectivdy.  The 
average  bit  error  probability  P,  is  obtained  by  taking  the  expected  value  of  P,  (b)  in  (18) 
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over  all  bit  patterns  b.  The  minimuin  bit  error  probability  is  (Stained  by  optimizing  over 
the  threshold  a.  In  summary,  given  a,  we  calculate  the  following  expectations; 


Pe  =E  {Peib)}  (20) 

b 

P,,^=  imn£  {/>.(*)>  =nim  /K*) 

where 


2ic  2% 


2n  2n 


+  J  •  •  •  J  •  •  •  d^MI2  } 

n  Af  0  V^O/ 


0  0 


(21) 


m.  NUMERICAL  RESULTS 


In  this  section  we  present  the  numerical  results  for  a)  hit  error  probability  versus 
signal-to-noise  ratio  Z  =  /fPjTINo  as  a  function  of  the  normalized  chaimd  spacing 
(normalized  to  the  bit  rate) 

where  4/*=  od^/  Ink  is  the  equal  channel  spacing  in  CiHz,  and  b)  power  penalty  versus 
normalized  channel  spacing  /  as  a  function  of  the  FP  filter  parameter  cT  =  {fd  Jp^BT 
where  B  is  the  filter  full  width  at  half  maximum  bandwidth(FWHM) 

Figures  4-5  show  the  minimum  bit  error  probability  versus  signal>to>noise  ratio 
(  Z )  for  cT  =  5  and  10,  respectivdy,  and  that  of  a  single  channd  (SC)  operation  without 
filtering  or  with  filtering  but  without  ISI.  In  Fig.  4  we  observe  that  a  large  degradation 
occurs  due  to  ISI  for  cT  =  5  which  represents  a  narrowband  filter.  As  the  FP  filter 
bandwidth  is  made  larger  as  in  Fig.  5  with  cT  =  10,  the  ISI  is  reduced  but  the  ACI 
increases. 

In  our  model,  we  are  constrained  to  Af  -  4  for  the  case  under  consideration.  We  set 
a  FP  fiber  with  fi’ee  spectral  range,  FWHM  ^=121 .6  GHz,  FSR  =  3800  GHz,  finess  F  = 
FSm  =  312.6,  and  c  =  38  4  GHz,  1/7'=  2  56  Gb/s,  then  cr=  15 
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Probability  of  Btt  Error 


signal-to-noise  ratio  as  a  function  of  normalized 


Figure  4;  Probability  of  bit  error  versus 
channel  spacing  for  cT-S 


\’i 


Probability  of  Bit  Error 


Figure  5;  Probability  of  bit  error  versus  signal-to-noise  ratio  as  a  function  of  normalized 
channel  spacing  for  cT=\0 
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(i  e  channel  spacing  is  12  times  the  bit  rate  or  30  4  GHz)  In  this  model  the  farthest 
adajcent  channel  for  A/  =  4  is  twice  the  channel  spacing  which  is  60  8  GHz.  This  verifies 
the  assumption  \f-fo  \  <  FSR/20n  =  60  5  GHz,  where  is  the  FP  filter  center  frequency 
This  result  agrees  well  with  that  in  [Ref  6,  Figs.  6,9,  MIF  =  0.4,  a  =  0.2].  Thus  we 
incorporate  the  degradation  caused  by  the  four  nearest  adjacent  channels.  We  observe 
that  for  bit  rates  of  1  Gb/s  or  less,  our  model  is  valid  up  to  M  =  10,  and  very  little 
difference  is  observed  between  M  =  4  and  A/  =  10  Also,  we  observe  that  there  is  little 
difference  between  M=2  and  A^  =  4  when  /  >  =10  for  bit  rates  up  to  3  Gb/s  In  all  results 
we  set  Lg  =  2  and  Z.  =  0 

Figure  6  also  shows  the  power  penalty  for  a  dense  WDM  system  relative  to  a  single 
channel  operation  at  the  minimum  bit  error  probability  of  10  '’  This  is  the  required 
additional  signal  power  (dBW)  for  the  dense  WDM  system  to  be  able  to  operate  at  the 
lO  ”  bit  error  probability  achieved  in  the  single  channel  system  with  a  SNR=12dB  The 
dense  WDM  system  is  /S7-limited  at  2  2  dB,  1  dB,  0.5  dB,  and  0  4  dB  in  power  penalty 
for  cT  =  5,  10,  15,  and  20,  respectively.  It  is  seen  that  for  a  2.3  dB  power  penalty,  the 
normalized  channel  spacing  can  be  as  close  as  /  =  6  (i.e.,  a  channel  spacing  of  six  times  the 
bit  rate)  for  cT=  5.  If  the  power  penalty  criterion  is  I  dB,  the  normalized  channel  spacing 
is  7=12  for  cT=  10,  15,  20.  We  remark  that  although  the  exact  transfer  functions  of  the 
FP  filter  is  used  in  [Ref  6],  a  number  of  approximation  have  been  made  to  obtain 
numerical  results.  The  approximations  are  1)  the  JSI  is  obtained  by  modeling  FP  filter  as 
a  single-pole  RC  filter  [Ref  6,  Eqs.  (4)  and  (36)],  2)  approximating  the  finite  integration 
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Power  Penalty  (dB) 


with  an  infinite  integration  in  the  calculation  of  ACI  [Ref.  6,  Eq.  (IS)],  and  3)  the  beat 
interference  is  ignored.  On  the  other  hand,  the  ISI  and  ACI  in  our  investigation  are 
obtained  by  modeling  the  FP  filter  as  a  single-pole  RC  filter,  using  finite  integration  and 
including  the  beat  interference.  Since  the  results  in  our  investigation  and  in  [Ref  6]  agree 
well,  we  conclude  that  approximations  are  quite  valid.  We  also  note  that  our  resuits  also 
agree  well  with  the  simulation  carried  out  in  [Ref  1,  Fig.  17]. 

The  above  numerical  results  shown  in  Figs.  4-6  are  obtained  with  an  optimized 
threshold  setting.  Figure  7  shows  the  power  penalty  for  fixed  threshold  a  =  £PTI1  which 
is  the  same  optimum  threshold  for  single  chaimel  operation  (midpoint  between  the 
received  power  for  bit  zero  and  bit  one)  It  is  seen  that  the  performance  of  a  dense  WDM 
system  is  quite  sensitive  to  a  for  a  narrow  band  filter.  An  additional  1.8  dB  is  observed 
for  c7  =  5  for/>  8,  and  0.5  dB  forcJ  =  10for/>  12.  Negligible  degradation  is  observed 
for  cT=  15,  20  for  7  >  16. 

Figures  8-9  show  the  power  penalty  versus  normalized  channel  spacing  as  a  function 
of  FP  filter  parameter  cT  for  the  worst-case  analysis  with  optimal  threshold  and  fixed 
threshold,  respectively.  The  worst-case  bit  pattern  is  fixed  to  produce  the  minimum  A', 
and  maximum  where  X■^  and  X^  are  the  values  to  A"  in  Appendix-A  equation  (9)  with 
=  1  and  Aqo  =  0  respectively 
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Figure  7:  Power  penalty  versus  normalized  channel  spacing  as  a  function  of  Fabry-Perot 
filter  parameter  cT  with  a  fixed  threshold  a=0.5 
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gure  8:  Worst-case  power  penalty  versus  normalized  channel  spacing  as  a  function  of 
Fabry-Perot  filter  parameter  cT  with  optimal  threshold 
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Power  Penalty  (dB) 


Figure  9;  Worst-case  power  poudty  versus  normalized  channd  pacing  as  a  function  of 
Fabry-Perot  filter  parameter  cT  with  fixed  threshold  a’^.S 
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We  observe  that  the  power  penalty  for  the  worst-case  analysis  is  only  slightly  larger 
than  that  of  the  exact  analysis  for  7  >  10  shown  in  Fig  6  Similarly  the  power  penalty  for 
the  worst-case  analysis  with  fixed  threshold  is  only  slightly  larger  than  that  of  the  exact 
analysis  with  fixed  threshold  for  7  >  10  shown  is  Fig.  7.  The  reason  for  this  is  that  for 
large  channel  spacing  (7  >  10),  the  ACI  effect  is  small,  so  the  ACI  bit  pattern  has  a  small 
influence  on  the  power  penalty. 

Figure  10  shows  the  normalized  optimal  threshold  for  the  exact  analysis  shown  in 
Fig.  6.  It  is  observed  that  a  «0.4  for  7  >  10.  Note  that  the  normalized  optimal  threshold 
for  the  single  channel  operation  is  a  =  0.5. 
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optimal  Threshold 


IV.  CONCLUSIONS 


We  have  preseitted  a  simple  modd  fyr  the  analysis  of  dense  WDM  systems 
employing  an  external  OOK  modulator.  The  only  ai^oximadon  that  we  use  involves  the 
modding  of  the  Fabry-Po-ot  fiher  by  a  sin^pole  RC  fiber  assuming  the  equivalem 
lowpass  signal  is  bandlimited  to  the  frequency  lange  |/|  <  FS/l  /  20x  This  modd 
enables  us  to  obtain  a  closed  form  expression  for  the  tm  error  probability  which  previously 
can  only  be  obtained  via  numerical  analysis  [Ref  6].  For  FP  filter  with  an  FSR  around 
3800  GHz,  our  nuxlel  can  include  the  AC/  effects  of  four  adjacent  channds  for  bit  rates 
up  to  2.S  Gb/s.  Our  numerical  results  show  that  this  modd  agrees  wdl  with  that  in 
[Ref  6] 
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APPENDK  A 


DERIVATION  OF  FORMULA  FOR  DECISION  VARIABLES 

DESIRED  CHANNEL  CHANNEL  0 

ADJACENT  CHANNEL  k,k  =  -A//2, “1, 1, A//2;  A/ :  even 
BIT  IN  CHANNEL  0  IN  ith  TIME  INTERVAL  (iT,  (/  +  1)7) 
haj  €  { 0, 1 } ,  hofi  :  DETECTED  BIT  IN  (0,  7) 

BIT  IN  CHANNEL  k  IN  Ith  TIME  INTERVAL  {IT,  (/  +  1)7) 
bkj  e  { 0,  }  WHERE  j  =  imaginary  numbler 

FREQUENCY  SPACING  BETWEEN  CHANNEL  k  AND  CHANNEL  0. 

^  PHASE  OFFSET  BETWEEN  CHANNELS  k  AND  CHANNEL  0. 

UNIFORMLY  DISTRIBUTED  BETWEEN  (0,  2jc) 

DATA  SIGNAL  IN  CHANNEL  0 

M0=  i  bojpTit-iT) 

i=-Lo 
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DATA  SIGNAL  IN  CHANNEL  k 

*t(0=  i  btjei^^piit-lT} 

L^,L  INTEGERS 
0<f<7' 

pTiJ)  “"10,  otherwise 


_  /a  -Twt  L  *T</<(i+l)r 

Pryt  -  !/)-{  0,  otherwise 

THE  RECEIVED  SIGNAL  AT  THE  INPUT  OF  THE  FP  FILTER  OF 
CHANNEL  0  IS 

r(/)  =  y?M/>+ T  JPbM 

k=-MI2 

M 

P  RECEIVED  OPTICAL  POWER 
THE  OUTPUT  OF  THE  FP  FILTER  IS 

ro(/)  =  f  hit-x)r(x)dx 

-00 

WHERE  MO  IS  THE  EQUIVALENT  LOWPASS  IMPULSE  RESPONSE  OF 
THE  FP  FILTER  OF  CHANNEL  0 

roit)=JP  j  h(t  -  x)bo(x)clx  +  JP  ^  J  h{t-x)bk{x)dx 


=  JPbo,o  f  h{t-x)pTix)dx 


k=-MI2Ac 

heO 


+Jp  sl  6o,/  /  hit-x)pTix-iJ)dx 

;=-Z,0  -OC 
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^  i  bk,i  f  h{t-x)e^^^^pjix-lT)dx 

k^-mt^L  4o 
M 

SINCE  THE  DETECTION  INTERVAL  IS  0  <  /  <  T,  WE  ONLY  NEED  TO 
EVALUATE  s(t)  =  rjit),  0<t<T 

s{t)  =  JPbQfi  /  /r(/-  t>/t  +  /?  2l  ^0,/  /  >*(/  -  x)dx 

0  /=-Lo  /T- 

(I- 

bkj 
k*0 

+^*,0  /  /rCZ-ty^^VT} 

0 

0<t<T 

S^(t)  DESIRED  SIGNAL 

S^/t)  INTERSYMBOL  INTERFERENCE 

S4cX0  ADJACENT  CHANNEL  INTERFERENCE 
FPITLTER 

LOWPASS  EQUIVALENT  TRANSFER  FUNCTION 


+Jp  1  {[^ 


*=^W/2 


h  )T 

h(t-x)e^^'dx] 
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31 


00 


J 

-00 


h{t-z)pr(T)  =  0 


00 

J  h(t-x)priT)dT 


00 

J  h{t-T)pT(T)dx 


-oc 


= J  /?(/  -  z)dx 

0 


=J  h{t-x)dx 

0 
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(3) 


Sait)  =  JPbo^o  J  h{t  -  x)dx  0<t<T 

0 

SUBSTITUTING  (2)  INTO  (3)  WE  OBTAIN 

SB(t)  =  JPbo^Q  /  ce-^*-^^dx  =  JPbo,oce-^'  /  e^iH 
0  0 

=  JPbo.oce-^'y^\o  =  1) 

=  JPbo,o(\  -  ^-^0  0<t<T 

a)  BIT  1  5b,i(/)  =  JP(l  -  0<t<T 

b) BIT0  5B,o(/)  =0  0</<7' 

^  ^isi  • 


J  h(t -  T)pr(‘c  ~  =  0 

-00 
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» 

in 

1 

/  >  (/  + 

l)T 

B 

iT  (1+1)7- 

V _ 

“  / 

> 

J  h{t-x)pTix-iTydx 

-00 

=J  h{t-x)dx 

iT 


J  h(t-x)prix-iT)cix 

-00 

(1+1)7- 

=  J  h{t-x)dx 

iT 


Sisi{t)=jp  t  boj  /  h{t-x)dx  0<t<T  (5) 

i=-Lo  iT 


SUBSTITUTING  (2)  INTO  (5)  WE  OBTAIN 


^  -1 

sisi(t)  =  jp  2  boj  J  ce-^'-'^dx 

l=-U  IT 


=jp  i 

i--Lo 


sisi(.t)  =  jp  i 

i=-Lo 
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sis,(l)  =  JPe-"  ^ 

i=-Lo 

WORST-CASE  ISI  boj  =  b-  ,  Lq 

=  V?  b- 

3.  s^cXt) : 


SUBSTITUTING  (2)  INTO  (7)  WE  OBTAIN 


T_1 

{[i  btj  J 

k=-M/2  l^-L  IT 

k*Q 

+*t,o  J  ce-^'-'V“‘VT} 

0 

A//2  —1  (^'*■1)7' 

=  JP  Z  {[2  bijce-“  j 

*=-M2  /=-!  IT 

k*0 


+bt,oce-<^  J 


0 


=  jp  f  {[  ^  btjce 

k=-M/2  l=-L 
k*0 


+bk,oce 


C+J&k 


-1 


SAcO)  =  JPce-'  ^{  [ 


o</<r 


WORSE-CASE bkj  =  bk,o  =  b,L  =  oo 


g(c+j<ak)lTy^ 


(8) 
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*=-M/2  c(l+/-r) 
iM) 


SAcAf)=JPf>  ^  0</<7’  (8.) 

'  fcr-A//2  1+y-r 

k*0 

OUTPUT  OF  PHOTODETECTOR  ;  ^  1^(0 1  ^ 

T  2 

SIGNAL  OUTPUT  OF  INTEGRATOR  :  |5(/)|  dt 

0 

I  j(0l  ^  +  SlSliO  +  5y<C/(0l  ^ 

=  l•sfi(0l^  +  |5/s7(0P+I^><c/(0P +2/?^{j5(0*y/s7W} 
+2Re{s*B(i)sACi(0}  +  2/?e{5j5y(0^yiC/(0} 

WORSE^ASE 


\)SB{t)  =  JPbo,o(\-e-^^)  0<t<T 

l55(0l  ^  =  Pbioi^  -  =  Pblo{\  -  2e-^'  +  ^‘^0  0<t<T 

J  \sB{t)\^dt  =  0  -  2P^o,o  J e-^^dt+Pblo  ] 

0  ’00 
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=  PTbl^  -  lPbla^e^'\l + 


=  Pni,  +  -  1)  -  - 1) 

=  PTbl<A\  1)] 

=  PTbU^-jji\  -6-^0  + 2^1  -e-'^^)] 


2)  s“^)  =  JPb-e-^'  0<t<T 

kSWl' o<'<^ 


J  b^Wl^  t=Pbl]  e-^dt  =  Pbl:^e-^\l 
0  0 


=  -^(e-^^-l)  =  f(l-e-^0 


n^/\  -^-2cT^ 

2cT 


3)«;ja(o=v?6  T  o<<<r 

y|C/v  /  1+y-^ 

ib^ 


38 


ACi\/  *=-A//2/f=-A//2  (i+y-rXi-y-r) 


0<t<T 


M  h*0 


!  ■ '■'‘I  ’  X  iJwfcK  1 


M  h*0 


=/>iiP  T  f 


k=^/2h=^l2 
ifc^  feto 


M/2 


=  P\b\^  1 


f 


ic=-M/2  h=-Mi2 
lc*0  h*Q 


[ejQ-^^)cT^  1 


_  PT\b\^  ^  ^ 


€»&  «»J, 

eA-/-^)c7-., 


*=-M2  /i=^/2  y(  l+y^X  i-y T^X^r-^) 

fctO  /ntO 


©it  =  ^,  /  =  INTEOTR  >  0 


THEN  J  =  I 

0 


r,  c8k=o* 


^  —in 

i  *=-M/2  l-K-r)^ 

M 
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*)2Re{sB(l)s'a,(t))  =2Re{JP  bo.o(i -e-‘')JPb-e-‘'} 

=  2Pbo.ob-(e^'-e-^')  0<t<T 


]  2Re  {SBit)shi(t)  =  2Pbo.ob-[^\  J  -  ^1 J } 


-cr_i  , 


=  2/>^0,0^-{^  +  ^ 


PTbo.ob- 

cf 


(\  -26-^^  +  e-^^) 


5) 


2^e{iJ(0s^a(/)}  =  2Re{JPboA^  -e^)JPb  /f 


1+/-^ 
k*Q 


—  2Pbo^oRe{b 


2-  -  ®*  i 


k=-M{2  1+y-f 

htO 


] 2Re{s'B(.t)SAcM}dt  =  2Pbo.oRe{b  T 


*=-M/2  i+y- 
k*0 


CO  It  =  I  =  INTEGER>0 


0  0 
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h  -(c-j^,)  c  c 

J  2Re{sl(.t)SAcmdt  =  2Pbofilie{h  f  (^X-4^)> 

n  *s-A472  1+;—  ‘-V-r 


2fr&o.o^(^} 

cT 


(,_e-cr)  T  ^ 

I-K-r)^ 


6)  2Re{s',si(t)sAa(t))  =  2Re{JPb-e-«JPb  T  ^ 


fc=^/2  l+/~ 
M 


=2Pb-mb  f 

k=-M2  1+7-r 
ihtO 


J  2/ee{i;fl('W/)>-  2Pb-Re{b)  Jf  ^  i „-,2x,^:a, 


1  -e' 


T  -irT 

i=-^/2  H-r)^ 
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/=  INTEGER  >0 


SUMMARY  OO 


C0i  = 


iKkI 

T 


x=  if] 

0 


=  ^{/>riJ,o[i-^i +  2^1 


+pnb\ 


f 

M 


1-K?)^ 


+«^i 


-2e-^^+e-^^) 


IPTboflRejh) 

cT 


(l-e-'O 


f 

M 


2PTb~Re{b) 


cT 


f 


1 


k=:~MI2 


} 


+*-IS7<*  -  +  *o.o*-[^l  -  2e-"'’ + e-^0] 
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WHERE  {0,e^},/&{6}  e  {0,cos4},|i|^  e  {0, 1} 


EXACT  ANALYSIS 

X)\  M/)l'*  =  PraJ,otl  -^1  -e-*’')+5?(l  - 
0 

2)  sis(l)=  JPe-^  i  boAe^‘*'^’'-e“^ 

i=-Lo 

I^kKOI  ^  =  Pe-^'[  ^  boAe^'*'^’’  - 

i=-Lo 

]  U/sfWl  V/=  -e-^%  i  boAe^'*'^^ 

0  i=-Lo 


2)sacM  =  f 


k^M/2  c(l+y-r)  l=-L 
M 


0<t<T 

0<t<T 

^)(M)cr 


0<t<T 


sacM  =  JPe-<'  f  i 

k=-Mr2  h-L  l+/-r 
M 

+  _  e-") 

k=-M/2  1+y-r 
M 

0<t<T 


\sAcm\^  =  Pe-^'\  t  t 

'  '  k=-MI2l=-L  l+y-r 

M 


+/>  T  T  ".f’"..  (g»“«'-e^Xg'^‘'-g^') 

ii=-A//2  «=-A<72  (i+y-rXi-y-r) 

ffn«0 


+2Pfe{  f  f  i 

\h=-M/2  m=-M2  /=-!  (i+y-rXi-y*^) 

ibtO  ifn*0 


EVALUATE  /  |5,4C/(0P^*  • 


w 

^)]\sAcM\\dt=^\-e-^^)\  T  7^ 

j  I  /II  acT'  k^-Mni=-L  i+;-r 


(^d+y^XZ+Dcr  _  ^(i+y?)/cr)|  2 


SPECIAL  CASE  (0*  =  ^ 


/=INTEGER>0 


'  k=-M/2  h-L  »+yV 


vt]  WcMWdt  =  p  f 

J  I  >1WV/I2  m*-A//2  (i+j-rXH-r)  0 


ktO 


T  T  T 

_f  e-^‘-i^‘^dl-]  e-<‘*^-^dt+l e'^'dt } 
0  0  0 


_o  W  ^  «>*=«», 

—  /  i  i  .  tt-  I  J(»t-«m)7*  , 

ifc=;:M2  m=^/2  iH-rXH—) 


ib^  ifR^ 


.  e-<c»/>„)r.i  _i_  __  e~2cr\| 

C';/C()fc  C+TWrn  2c'‘ 


PT  i  iCT ,  (S}k=<am 

—  "7r  ^  "kwi  “w.  I 

k=-Mi2m=-Mi2  (1+i-rKi-;— )  ^  ^  ^-1 


k*0  m*0 


Ac  C  A 


^-(i-;^)cr.|  ,  _»-2cT\\ 

^  i^jr  2^  ^ 


2nkl 

SPECIAL  CASE  Ojfc  =  ^ 


/=  INTEGER  >0 


(Ow  = 
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0 


'  I,  (!)*«>, 


gJ&kT  _  g-MkT  _  I 


=  PTf  f  -M-^±  W  ^  ^>.o^:.o  re-^-l 

'as-^/2  k=-MI2  m=-M2 

ybto  Aseo  arK) 


.  1 


2 


+  i(l  -e-^^)]) 


C)]  \sAa(t)\ldt  =2PRe{  T  ^  ^ 

••  ^  ^  '  -J  ^  t_  t  _  ft  t  ft  /  I  .  .'_*  W 1  ^ 


fe-M2  m=-A//2  /=-!  (i+;-rKi-;^) 
jbtO  m?tO 


:Re{ 


M  m*0 


0 

0 

**./6:.o 

:i+y?Xi-;^) 

-in 

SPECIAL  CASE  (Ok 


/=INTEGER>0 
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M  nnK) 


-kl -e-^^)]] 


]\sAat)Vdt=^\-e-^H  f  i 

J  I  iici  /I  2cn  '  k=^ni=-L  i+;-r 


*=-A//2«=-A//2  (l+;-rXi->-r) 
ib^  m^ 


..cf,  <iiie~<am  g-(H 

I  I  ^ 

fUlc  c  ^‘~\  ^  t~J’T 

A-h^) 


+2^e{  T  T  ^ 

ik=-A</2  m=-A//2  /=-!  ( 1+y-rX  »-;-r) 
ifc^  m*0 


AlcT\t  ***'' 


c 


-id-e-^O]} 
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SPECIAL  CASE 


®  k*0 


^  '*=^/2  1-K^)^  /k=-M/2  m=-A//2  (1+/^X1-/ 

jbtO 


2ioit/  V 

cr  ^ 


rlZlzi.  +  i^Id.  +  In  -  I 

l-u.iaw  ^  i+,i25i  ^  2'^* 

^  J  rT  cf 


-2cr>i 


cT  -  cT 


+2^/te{  f  f  i 

lc=-M/2m=-mi^-L  (1+/^X1-;V> 
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APPENDIX  B 


MATLAB  COMPITTER  PROGRAMS 

exact  analysis  with  optimal  threshold  signal  out  of  the  integrator,  the  signal  X  cmttains  desired 
signal  and  ACI  &  ISI  and  postdetect  noise  .The  formula  for  the  BER  is ; 

M  is  number  of  adjacem  channels  Since  we  assume  upper  and  lower  channels  are  synchronized 
individually,  i.e.  for  our  model  the  power  of  (l/2ic)  is  fixed  to  2,  also  we  only  have  two 
integrations  and  two  arguments. 

M=4,  k=[.M/2>l  l:M/2], 
m=k; 

%  produce  the  controlled  matrix  b  to  control  64  difierent  bit  patterns 
ml=[  zeros(l,32)  ones(l,32)  ]; 
m2=[  zeros(  1,16)  ones(  1,16)  zeros(  1,16)  ones(  1,16)]; 
m3=[]; 
for  i=l  4 

m3=[m3  [zeros(l,8)ones(l,8)]]; 
end 
m4=[], 
for  i*l;8 

m4*[m4  [zeros(l,4)  ones(l,4)]]; 
end 

m5*zeros(l,64); 
m6=[]; 
for  i=l:16 

m6=[m6  [zeros(l,2)  ones(l,2)]]; 
end 
m7=[], 
fori=l  32 
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m7-{m7  zeros(l,l)oiies(l,l)]. 
end 

b-fml  ,in2;m3;iii4;m5;ni6;m7]'; 

%  signal  to  noise  ratio  range  in  dB 
RPTN  DB»[10:0  5:25];  •/•RPTN_dB 
ppp*  1 0  /XO  1  ♦RPTN_DB); 
lenl=length(ppp); 

%  solalph  function  is  provided  by  Professor  Randy  L.  Borchardt 
n*10; 

[bpx,wfx]*gnile(n);  %bpx*4)py  ,wfx=wfy 

%  single  channel 
BER0=O  5*erfc(ppp/8^.5). 

pp=[],  threslr»[];  %thresh  is  not  nomudized  threshold 

for  CT-[S  10  IS  20]  %  cT  is  Fabry-Perot  filter  parameter 
ifCT«=5 
qqq«(2:12  17  20]; 
elseifCT=10 
qqq=[3:12  17  20]; 
elseifCT=15 
qqq=[4:12  17  20], 
elseifCT='20 
qqq=[5:12  17  20]; 
end 

pe=[];  threshl=[]; 
for  I=qqq 

BER=[];  thresh2=[]; 
for  RPTN=ppp 

ap^’^Unspat^O,  1 82, 1 1 );  Voapproximated  thresholds 

}b=0;  %  fim  finv  loops  to  find  out  the  threshold  which  nuke  the 

%BERmininnin  don't  care  about  the  scale 

fori=l:64  x3=x3+... 

solalph('xx00',0,2*pi,2,bpx,w6c,0,2*pi,2,bpx,wfic,ap,b,CT,i,I,k,m,RPTN)+. . . 
solalph('xxir,0,2*pi,2,bpx,wfx,0,2*pi,2,bpx,wfic,iq),b,CT,i,I,k,m,RPTN); 
end 

[val,ind]>min(x3); 


ief^ap(ind)-16  6, 

if  lef^O  ,lef<K};  end  %  to  aviod  the  threshold  go  beyond  the  negative  side 

i^linspace(l^ap(ind)+16.6,l  1), 

%  two  more  time  to  find  alpha 
x3=0; 


fori=l;64  x3=x3+... 

solalph('xx00'.0,2*pi,2,bpx,w6t,0,2*pi,2.bpx,wfit,ap,b,CT,i,I,k,m,RPTN)+ 

soialph('xxll’,0,2*pi,2.bpx,wfit,0,2*pi,2.bpx,w6c,ap,b,CT,i,I.k,m,RPTN). 

end 

[val,ind]=min(x3); 

lef=ap(ind)-3.1, 

iflef<0,lefN);  end 

q}=linspace(lef,ap(ind)+3 .1,11); 

%  three  more  time  to  find  alpha 

x3=0; 


for  i*l;64 
x3=x3+.. 

solalph('xx00',0,2*pi,2,bpx,wfx,0,2*pi,2,bpx,wfit,ap,b,CT,i,I,k,m,RPTN)+ 

solalph(’xxir,0,2*pi,2,bpx,w6t,0,2*pi,2,bpx,w6i,ap,b,CT,i,I,k,m,RPTN); 

end 

[val,ind]=min(x3); 

lef^ap(ind)-0.57; 

if  lef^O  ,lef=0;  end 

ap=linspaceOef,ap(ind)+0. 57,8); 

%  four  more  time  to  find  alpha 
x3=0; 


for  i=l;64 
x3=x3+... 

solalph('xx00’,0,2*pi,2,bpx,w6c,0,2*pi,2,bpx,wfic,ap,b,CT,i,I,k,in,RPTN)+ 

solalph('xxll’,0,2*pi,2,bpx,wfx,0,2*pi,2,bpx,wfii,ap,b,CT,i,I,k,m,RPTN); 

end 


[val,ind]»nnin(x3); 

ap=ap(ind); 
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%  after  find  optimal  alpha  use  double  integration 
qq3»0; 

fori=l:64 

qq3=qq3+(  dbgquadni('xx00',0,2*pi,2,bpx,wfic,0,2*pi,2,. . 

bpx,A^ap,b,CT.i.I,k,m,RPTN)+. 
dbgquadniCxxl  r,0,2*pi,2,bpx,wfi£,0,2*pi,2,. . 
bpx,w6c,ap,b,CT,i,l,k,m,RPTO)  )/( 1024*pi^2); 
end 

BER^'CBER  qq3];  thresh2»[thresh2  ap]; 

if  qq3<l(yX-15)  %for  save  time  only  interesting  in  10^-15) 
ii=find(ppp*=RPTN); 

BER(ii+l:lenl)=5*l(KX-l  16)*ones(  l,length(ii+l;lenl)  ); 
thresh2(ii+l:lenl)=(*P+2)*ones(l,lei!®th(ii+l;lenl)  ); 
break 
end 


pe=[peJBER];  thresh l=(threshl;thresh2]; 
end 

pp=lpp;pe];  thresh*[thresh;thrcshl]; 
end 


tinic2=toc. 


%  worse  case  form  appendix  setting  optimal  threshold  equal  (x0+xiy2 
M*4, 

I=linspace(0,6, 101);  %  set  I  value  in  linspceinteger  to  find  minimum  I 

RPTN_DB=0:0.2;20,  %  x-axis  dB  range 

RPTN=  1 0/XO  1  *RPTN_DB);  %  change  to  ratio 

%  single  channel 

BER0=0  5  *erfc(RPTN/8'H)  5); 

%find  optimal  alpha 
PP=[]. 

forCT=[5  10  15  20] 
x3-0; 

fork=l:0  5*M 

x3=x3+2 /(l+(2*pi*k*I/CT)/^2), 
end 

x0=(  1  -exp(-2*CT))/(2*CT)+x3*(  1  +2*(  1  -exp(-CT))/CT). 

X 1 = 1  -2  *(  1  -exp(-CT))/CT-K  1  -exp(-2  *CT))/(2*CT), 
[val,ind]=min(abs(xO-xl )), 

I=ceil(I(ind)), 

I=(I(ind)),  %  our  minimum  I  value 

%  set  different  I  and  to  find  the  BER 
ss=H-l:I+10, 

BER*[], 
for  I=ss 
x3=0; 

for  k=l:0.5*M 

x3=x3+2/(  1 +(2*pi*k*I/CT)^2), 
end 

x(H  1  -exp(-2*CT))/(2*CT)+x3  *(  1  +2*(  1  -exp(-CT))/CT); 

X 1 = 1  -2*(  1  -exp(.CT))/CT+(  1  -exp(-2*CT))/(2*CT); 

BER=[BER;0  5*erfc(  RPTN'»(xl-x0)/8'X).5  )  ]; 
end 

pp=[pp;BER]; 


semilogy(RPTN_DB,BERO'~',RPTN_DB,pp(;,;)) 
axis([10  19  lO'X-lS)  1]) 
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